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Abstract 



The Fornax galaxy cluster was observed with the High Energy Stereoscopic 
System (H.E.S.S.) for a total live time of 14.5 hours, searching for very-high- 
energy (VHE, E > 100 GeV) 7-rays from dark matter (DM) annihilation. No 
significant signal was found in searches for point-like and extended emissions. 
Using several models of the DM density distribution, upper limits on the DM 
velocity-weighted annihilation cross-section (av) as a function of the DM par- 
ticle mass are derived. Constraints are derived for different DM particle mod- 
els, such as those arising from Kaluza-Klein and supersymmetric models. Var- 
ious annihilation final states are considered. Possible enhancements of the DM 
annihilation 7-ray flux, due to DM substructures of the DM host halo, or from 
the Sommerfeld effect, are studied. Additional 7-ray contributions from inter- 
nal bremsstrahlung and inverse Compton radiation are also discussed. For 
a DM particle mass of 1 TeV, the exclusion limits at 95% of confidence level 
reach values of cm^s , depending on the DM particle 
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model and halo properties. Additional contribution from DM substructures 
can improve the upper limits on (av) by more than two orders of magnitude. 
At masses around 4.5 TeV, the enhancement by substructures and the Sommer- 
feld resonance effect results in a velocity-weighted annihilation cross-section 
upper limit at the level of ^(7-^^95% c.L. r^\Q-'2.6 cm'^s"^. 

Subject headings: Gamma-rays : observations - Galaxy Cluster, Dark Matter, 
Fornax galaxy cluster 



1. Introduction 



Galaxy clusters are the largest virialized objects observed in the Universe. Their 
main mass component is dark matter (DM), making up about 80% of their total mass 
budget, with the rer nainder pr ovided by intracluster gas and galaxies, at 15% and 5% 
respectively (see e.g IVoitI 120051) . The DM halo distribution within galaxy clusters ap- 
pears to be well reproduced by N-body numerical simulations for gravitational structure 
formation jColafrancesco et al\ bood: iRichtler et fl/l boosi: ISchuberth et al\ boid IVoitll2005 , 
and references therein). This may be in contrast to smaller systems like dwarf galax- 
ies. For instance, disagreements between theoretical predictions and actual estimates of 
the DM halo profile from observations have been found in lo w surface brightness galax- 
ies (|McGaugh and de Bloklll998l:lNavarrolll998l:lde Bloklboiol) . Although such discrepan- 
cies may vanish at galaxy cluster scale, the influence of baryon infall in the DM gravita- 
tional potential can still fl atten the DM dens ity distribution in the inner regions of galaxy 
clusters (see, for instance, lEl-Zant et fl/.ll200ll) . 



The pair annihilation of weakly interacting massive particles (WIMP) constituting 
the DM halo is predicted to be an important source of non-thermal particles, including 
a significant fraction a s photons covering a broad multiwavelen gth spectrum of emis- 
sion (see, for instance, iBergstromI l2000l: IColafrancesco et 120061) . Despite the fact that 
galaxy clusters are located at much further distances than the dwarf spheroidal galaxies 
around the Milky Way, the higher annihilation luminosity of clusters make them compa- 
rably good targets for indirect detection of dark matter. The flux of 7-rays from WIMP 
DM annihilation in cl usters of galaxies is possibly large e nough to be detected by cur- 
rent 7-ray telescopes (|Teltema et flZ.ll2009l:lPinzke et fl/.ll2009l). Also standard astrophysical 
scenarios have been proposed for 7-ray emission (see e.g. iBlasi et fl/.ll2007l, for a review), 
in particular, collisions of intergalactic cosmic rays and target nuclei from the intracluster 
medium. Despite these predictions, no signific ant 7-ray er riission has been obs erved in lo- 
cal clusters by H.E.S.S. jAharonian et fl/.ll2009allcl) , MAGIC jAleksic et fl/.lboiOah and Fermi- 
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LAT jAckermann et fl/iboiOallbl) collaborations. However, 7-rays of a different astrophys- 
ical emissior i processes have already been detected frq r n some central radio galaxies in 



cluste rs (e.g. lAharonian et all (|2006al) : lAcciari et all (|2008l) : lAleksic et all (|2010b[ ): lAbdo et al. 
(|2009|) V 



Following the absence of a signal, upper limits for a DM annihilation signal com- 
ing from galaxy clusters hay e been published by Fermi-LAI (|Ackermann et fl/.ll2010al) and 
MAGIC (jAleksic et alluOlOm collaborations. Strong constraints on the annihilation cross- 
section of DM from the Fornax galaxy cluster have been put by the Fermi-LAJ collab- 
oration for DM particles masses up to 1 TeV from 7-ray selected in the 100 MeV - 100 
GeV energy range. However many DM models show distinct features in the DM anni- 
hilation spectrum close to DM particle mass, such as monochromatic gamma-ray lines, 
sharp steps or cut-offs, as well as pronounced bumps. This could provide a clear distinc- 
tion between an annih ilation signal and a standard astrophysical signal (see, for instance, 
Bringmann et alluOllh ). These features are often referred as smoking-gun signatures. Such 
models can only be tested by satellite telescopes for DM particle masses up to a few hun- 
dreds of GeV. lACTs observation can provide well-complementary searches for such fea- 
tures at DM particle masses higher than a few hundreds of GeV . 

This paper reports on the observation in VHE 7 rays of the Fornax galaxy cluster 
(AGO S373) with the High Energy Stereoscopic System (H.E.S.S.). Interdependent con- 
straints on several DM properties are derived from the data, such as the DM particle 
mass and annihilation cross-section. Different models of the DM density distribution of 
the cluster halo are studied. The paper is structured as follows. In Section 2 the Fornax 
galaxy cluster is described. The choice of Fornax for a DM analysis is motivated, based 
on the DM content and distribution inside the cluster. Section 3 presents the data analysis 
and results. Upper limits on the 7-ray flux for both standard astrophysical sources and 
DM annihilation are extracted in Section 4. Exclusion limits on the DM annihilation cross- 
section versus the particle mass are given in Section 5. Several DM particle candidates are 
considered, with particular emphasis on possible particle physics and astrophysical en- 
hancements to the 7-ray annihilation flux. 



2. Target selection and dark matter content 



The Fornax (distance = 19 Mpc, Tonrv et al. II2OO1I), Coma (distance = 99 Mpc lReiprich and Bohringer 

2002) and Virgo (distance = 17 Mpc, Mei et flZ.ll2007 ) galaxy clusters are in prin ciple promis- 



ing ta rgets for indirect dark matter searches through 7-rays, as was shown by lleltema et al. 



(|2009l) . The radio galaxy M 87 at the center of Virgo provides a strong astrophysical 7-ray 
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signal (|Aharonian et fl?.ll2006ah , showing flux variabilities from daily to yearly timescales 
that exclude the bulk of the signal to be of a DM origin. Since a DM 7-ray signal would 
be hard to disentangle from this dominant standard astrophysical signal, Virgo is not a 
prime target for DM searches, even though a DM signal may be hidden by the dominant 
7-ray signal from standard astrophysical sources. 

Moreover, galaxy clusters are expected to harbor a significant population of relativis- 
tic cosmic-ray protons originating from diff erent sources, such a s large-sc ale shocks asso - 
ciated with acc retion and merg er processes ( Colafrancesco and Blasill 19981: Ryu et al. 



or supemovae jVolk et flZ.lll996|) and AGN activity (jHinton et al.\ 



2007|) . The 7-ray emission 



arising from pion decays produced by the interaction of these cosmic-ray protons with the 
intracluster gas may be a pot ential astrophysical background to the DM-tnduced 7-ray 
signal. In the case of Coma, iTeltema et al\ (|2009l) showed that such astrophysical back- 
ground is expected to be higher than the DM annihilation signal^ On the other hand, 
the same study ranked Fornax as the most luminous cluster in DM-induced 7-ray emis- 
sion among a sample of 106 clusters from the HIFLUGCS catalog (jReiprich and Bohringer 
20021) . The DM- to-cosmic-ra y 7-ray flux ratio of Fornax was predicted to be l arger than 
100 in the GeV energy range (jTeltema et al. II2OO9I) . A recent independent study by lPinzke et al 
(|201l|) has also predicted Fornax to be among the brightest DM galaxy clusters with a 
favorably-low cosmic-ray induced signal. Although the central galaxy of the Fornax clus- 
ter, NGC 1399, is a radio galaxy and could in principle emit 7-rays , the super-rnassiy e 
black hole at the center of this galaxy have been shown to be passive (jPedaletti et al. I l201l|) . 
Indeed recent observations of several clusters with the Fermi-LAT detector have shown 



no 7-ray signal (jAckermann et flZ.ll2010bl) , and the most stringent limits on dark matter 
annihilation were derived from the Fornax observations (jAckermann et flZ.ll2010al) . 



The center of Fornax galaxy cluster is located at RAQ2000.0) = 03^38"^29S3 and 
Dec(J2000.0) = —35° 27' 00"7 in the Southern Hemisphere. For ground-based Cherenkov 
telescopes like H.E.S.S. (cf. Section [S]), low zenith angle observations are required to guar- 
antee the lowest possible energy threshold and the maximum sensitivity of the instru- 
ment. Given the location of H.E.S.S., this condition is best fulfilled for Fornax, compared 
to the Virgo and Coma clusters. Therefore, Fornax is the preferred galaxy cluster target 
for dark matter searches for the H.E.S.S. experiment. The properties of its dark matter 
halo are discussed in more details in the following section. 



Also the two brightest radio galaxies, NGC 4874 and NGC 4889, lying in the central region of Coma 
may be potential sources of a standard astrophysical 7-ray signal. 
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2.1. Dark matter in the Fornax galaxy cluster 



The energy-differential 7-ray flux from dark matter annihilations is given by the fol- 
lowing equation: 



d<l>^(An, 
dE^v 



{(TV) 



8n m^j^ dEj 



X /(An)AQ, 



(1) 



where (av) is the velocity-weighted annihilation cross-section, the mass of the DM 
particle and dNy/dEy the photon spectrum per annihilation. The factor 



(2) 



reflects the dark matter density distribution inside the observing angle AO. The annihi- 
lation luminosity scales with the squared dark matter density p^, which is conveniently 
parametrized as a function of the radial distance r from the center of the astrophysical ob- 
ject under consideration. This luminosity is integrated along the line of sight (LOS) and 
within an angular region AO, whose optimal value depends on the dark matter profile of 
the target and the angular resolution of the instrument. 

Numerical simulations of structure formation in the ACDM framework predict cuspy 
dark matter halos in gal axies and clusters of galaxies ( Navarro et fl/.lll996 : Fukushige and Making 
19971: iMoore et all Il998f ). A prominent parametr ization of such halos is the "Navarro- 
Frenk- White" (NFW) profile (|Navarro et flZ.lll997|) , characterizing halos by their scale ra- 
dius Vs at which the logarithmic slope is dkijO/dlnr = —2, and a characteristic density 
Ps = 4|0(rs). This profile was shown to be consistent with X-ray observations of the intr- 
acluster medium of galaxy clusters. The DM density profile is given by: 



Pnfw('") 



1 + ^ 



(3) 



Another prediction of ACDM N-body simulations is an abundance of halo substructures, 
as will be detailed in section IZ2l On the other hand, in scenarios where the baryon tnfall 
in the DM gravitational potential efficiently transfers energy to the inner part of the DM 
halo by dynamica l friction, a flatteni ng of the density cusp into a core-halo structure is 
predicte d (see e.g. lEl-Zant et al. 2001 ). These halos can be parametrized by the "Burkert 



profile" (|Burkertlll996l) : 



PB{r) 



Po^c 



(r + rc)(r2 + r^) 



(4) 
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Again, the dark matter density falls off as ^ outside the core radius Vc, but it ap- 
proaches a constant value po for J' ^ 0. In the following, dark matter halos of both types 
are considered. 

A commonly-used approach for the determination of the DM halo in galaxy clus- 
ter comes from X-ray measur ements of the gravitationally bound hot intracluster gas. 
From the HIFLUGCS catalog (jReiprich and B6hringerll2002h , the virial mass and radius 
of Fornax are found to be Mvir ^ 10^^ M© and Rvir ^ 1 Mpc (corresponding to about 
6° in angular diameter), respectively. Under the assumption of a NFW halo profile in 
ACDM cosmology, a rela tion between the virial mass and the concentration parameter 
c = Rvir/ '"s was found by lBuote et al\ (|2007l) . The halo parameters can thus be expressed 
in terms of ps and fg and are presented in Table [T] This model is hereafter referred as 
to R B02. A similar procedur e was applied in the Fermi-hAI DM analysis of galaxy clus- 
ters (jAckermann et alWlOlOm . 



A different approach is to use dynamical tracers of the gravitational potential of 
the cluster halo, such as stars, globular clusters or planetary nebulae. This method is 
limited by the observability of such tracers, but can yield less model-dependent and 
more robust modeling of the DM distribution. However, some uncertainty is introduced 
by the translation of the tracer's velocity dispersion measurement into a mass profile, 
which usually implies solvi ng the Jeans equations under some simplifying assumptions 
(jBinney and Tremaind 120081) . From velocity dispersion measurements on dwarf galaxies 
obser ved up to about 1.4 Mpc, a dynamical analysis of the Fornax cluster by lDrinkwater et al. 
(|200l|) constrained the cluster mass. The associated DM density profile, hereafter referred 
as to DWOl, can be well described by a NFW profile ( Richtler et al. 200^ with parameters 
given in Table [T] 



Richtler et all (|2008l) have analyzed the DM distribution in the inner regions of Fornax 



by using the globular clusters as dynamical tracers. This allowed an accurate DM mass 
profile measurement out to a radial distance of 80 kpc from the galactic cluster centre, cor- 
responding to an angular distance of ^ 0.25°. The resulting velocity dispersion measure- 
ments can be well fitted by a NFW DM halo profile with parameters given in Tabled] This 
density profile (hereafter referred as to RS08) determination is in good agreement wi th the 
determination inferred from ROSAT-HRI X-ray measurements (IPaolillo et fl/.ll2002h. De- 
tailed analysis using subpopulations of globular clusters done in ISchuberth et al\ (|2010h 
showed that both a NFW and a Burkert DM halo profiles can equally well fit the globular 
cluster velocity dispersion measurements. Representative DM halo profiles using differ- 
ent sets of globular clust ers samples, hereafter referred as to SRIO a6 and SRIO aio, are 
extracted from Table 6 of ISchuberth et al\ (|2010l) . The parameters for both the NFW and 
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Burkert DM halo profiles are given in Table [H 

Using the dark matter halo parameters derived from the above-mentioned meth- 
ods, values of / were derived for different angular integration r adii. The point-spread - 
function of H.E.S.S. corresponds to an integration angle of ^ 0.1° dAharonian et flZ.ll2006b|) . 
and most often the smallest possible angle is used in the search for dark matter signals in 
order to suppress background events. However, since a sizable contribution to the 7-ray 
flux may also arise from dark matter subhalos located at larger radii (see Section |Z2)) , in- 
tegration angles of 0.5° and 1.0° were also considered. The choice of the tracer samples 
induces a spread in the values of the astrophysical factor J up to one order of r nagni tude 
for an integration ang le of 0.1°. Note that the measurements of iRichtler et al. (|2008h and 
Schuberth et all (|2010( ) trace the DM density distribution only up to 80 kpc from the center. 
In consequence the derived values of the virial mass and radius are significantly smaller 
than those derived from X-ray rn easurements on larger distance scales (see for instance 
figure 22 of ISchuberth et aUuOlw . Thus the DM density values may be underestimated 
for distances larger than about 100 kpc. On the other hand, it is well known that for 
an NIW profile about 90% of the DM annihilation signal comes from the volume within 
the scale radius r^. Therefore, even for NFW models with large virial radii such as RB02 
and DWOl, the main contribution to the annihilation signal comes from the region inside 
about 98 kpc and 220 kpc, respectively. 



2.2. Dark matter halo substructures 



Recent cosmological N-bod y simulations, such as Aquarius (|Springel et fl/.ll2008l) and 
Via Lactea ( Diemand et al. 20081) , have suggested the presence of dark matter substruc- 
tures in the form of self-bound overdensities within the main halo of galaxies. A quantifi- 
cation of the substructur e flux contributiori to the total 7-ray flux was computed from the 
Aquarius simulation by IPinzke et al\ (|2009l) using the NFW profile RB02 as the DM den- 
sity distribution of the smooth halc|^ The substructure enhancement over the smooth host 
halo contribution along the line of sight is defined as Bg^b ( AO) = 1 + £sub ( AO) / £sm ( AO), 
where >Csjn/sub(An) denotes the annihilation luminosity of the smooth host halo and the 



■^This halo is also well suited with respect to the others discussed in Sectior lZT] since substructures in 
the form of gravitationally bound dwarf galaxies to Fornax are observed up to about 1 Mpc. They are thus 
included within the virial radius predicted by the RB02 profile (Ryir — 1 Mpc). 
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/(AO) [lO^i GeV^cm- 










NFW profile 




Model 


rs [kpc] 


Ps [Mqpc 3] 


^max - 0.1° ^max - 0.5° ^max - 1-0° 


RB02 


98 


0.0058 


112.0 6.5 


1.7 


DWOl 


220 


0.0005 


6.2 0.5 


0.1 


RS08 


50 


0.0065 


24.0 1.2 


0.3 


SRIO aio 


34 


0.0088 


15.0 0.6 


0.1 


SRIO a6 


200 


0.00061 


7.0 0.5 


0.1 








Burkert profile 




Model 


Vc [kpc] 


Pc [M0PC-3] 


^max - 0.1° ^max - 0.5° ^max - 1-0° 


SRIO aio 


12 


0.0728 


15.0 0.6 


0.2 


SRIO a6 


94 


0.0031 


2.4 0.5 


0.1 



Table 1: Dark matter halo models for the Fornax galaxy cluster. The first three columns 
show the selected profiles discussed in Section IZTI with their respective NFW or Burkert 
halo parameters. The last three columns show the astrophysical factor /, calculated for 
three different integration radii. 



additional contribution from substructures, respectively. The former is defined by: 



^sm/sub(AO) = AO X /sin/sub(An) 



dO 



AO 



1.0. s. 



dl X p. 



sm/sub 



(5) 



where Psm/sub is the DM density distribution of the smooth halo and substructures, re- 
spectively. In order to perform the LOS integration over the subhalo contr i butio n, an 
effective substruc ture density p^uh is parametrized following ISpringel et al\ (|2008l ) and 
Pinzke et al\ (|2009h as: 



Psub(0 



A(r) 0.8C/:sm(R. 
47rr2R„,v 



?(r) 



where 



and 



A(r) = 0.8 - 0.252 ln(r/Rvir) 



Bir) = 1.315 - 0.8(r/R, 



-0.315 



(6) 



(7) 



(8) 



£sin(Rvir) is the smooth halo luminosity within the virial radius Rvir- The normaliza- 
tion is given by C = (Mmin/Miijn)'^'^^^, where Mmin = lO^M© is the minimum sub- 
structure mass resolved in the simulation and is the intrinsic limiting mass of sub- 
structures, or free-streaming mass. A conventional value for this quantity is Myo^ = 
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10~^Mq (jDiemand et al\\2006i) , although a rather broad range of v alues, dowri to M ]\m = 
10~^^Mq, is possible for different models of particle dark matter (lBringmannll2009h. As- 
suming a specific DM model, a constraint on M\[^ was derived by IPinzke et al\ (|2009h 
using EGRET 7-ray upper limits on the Virgo cluster and a lower bound was placed at 
^lim = 5 X 10~^ Mq. Nevertheless, the effect of a smaller limiting mass is also investi- 
gated in this work. 

Figure [1] shows the substructure enhancement Bgub over the smooth halo as func- 
tion of the opening integration angle. At the distance of Fornax, integration regions 
larger than ^ 0.2° correspond to more than 65 kpc. Beyond these distances the sub- 
structure enhancement exceeds a factor 10. This justifies extended analyses using inte- 
gration angles of 0.5° and 1.0°. Two values of the limiting mass of substructures are 
used: M\i^ = 1O~^M0 and Miim = 5 x 1O~'^M0, inducing a high and a medium value 
of the enhancement, respectively. The values of Bgub for the opening angles of 0.1°, 0.5° 
and 1.0° and for both values of Mu m are given in Table |2l These values are larger than 



those derived in lAckermann et all (l2010al). In their stu dv the substructure enhancement 
is calculated from the Via Lactea ( Diemand et al. 2008|) simulation , where a differe nt con- 
centration mass relation is obtained. For a careful comparison see lPieri et all ([201 11) . 



^max 


0.1° 


0.5° 


1.0° 




= lO-^M© 


4.5 


50.5 


120 




5 X lO-^M© 


1.5 


8.2 


18.3 



Table 2: Enhancement Bgui-, due to the halo substructure contribution to the DM flux, for 
different opening angles of integration ^max- The enhancement is calculated for two lim- 
iting masses of substructures and over the smooth DM halo RB02. 



3. Observations and data analysis 

The High Energy Stereoscopic System (H.E.S.S.) consists of four identical imaging at- 
mospheric Cherenkov telescopes. They are located in the Khomas Highland of Namibia 
(23°16'18" South, 16°30'00" East) at an altitude of 1800 m above sea level. The H.E.S.S. ar- 
ray was designed to observe VHE 7-rays through the Cherenkov light emitted by charged 
particles in the electromagnetic showers initiated by these 7-rays when entering the at- 
mosphere. Each telescope has an optical reflector co nsisting of 382 rouri d facets of 60 cm 
diameter each, yielding a total mirror area of 107 m^ jBemlohr et alhoo^ . The Cherenkov 



light is focused on cameras equipped with 960 photomultiplier tubes, each one subtend- 
ing a field of view of 0.16°. The total field of view is ^5° in diameter. A stereoscopic 
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reconstruction of the shower is applied to retrieve the direction and the energy of the 
primary 7-ray. 

Dedicated observations of the Fornax cluster, centered on NGC 1399, were conducted 



in fall 2005 (jPedaletti et flZ.ll2008h . They were carried out in ivobble mode (jAharonian et al 



2006bl) , i.e. with the target typically offset by 0.7° from the pointing direction, allowing 
simultaneous background estimation from t he same field of view. T he total data passing 
the standard H.E.S.S. data-quality selection (jAharonian et flZ.ll2006b|) yield an exposure of 
14.5 hrs live time with a mean zenith angle of 21°. 



m 



The data analysis was perfo rmed using an improved model analysis as described 



de Naurois and Rollandl (120091), with independent cross-checks performed with the 



Hillas-type analysis procedure (jAharonian et all l2006b[) . Both analyses give compatible 
results. Three different signal integration angles were used , 0.1°, 0.5° and 1°. The cosmic- 
ray background was estimated with the template model (jRowelll 120031) , employing the 
source region, but selecting only hadron-like events from image cut parameters. 

No significant excess was found above the background level in any of the integra- 
tion regions, as visible in Fig. |2]for an integration angle of 0.1°. An upper limit on the 
total number of observed 7-rays, , was c alculated at 95% confidence l evel (C.L.). 

The calculation followed the method described in lFeldman and CousinsI (jl998l) , using the 
number of 7-ray candidate events in the signal region Nqn ^rid the normalized number 
of 7-ray events in the background region Nqff- Since the normalization is performed 
with respect to the direction-dependent accepta i ice an d event rate, the background nor- 
malization factor for Nqff defined in lRowelll (j2003j ) is a; = 1. This is equivalent to the 
assumption that the uncertainty on the background determination is the same as for the 
signal, allowing a conservative estimate of the upper limits. This information is summa- 
rized in Table |3l 

A minimal 7-ray energy (Emin) is defined as the energy at which the acceptance for 
point-like observations reaches 20% of its maximum value, which gives 260 GeV for the 
observations of Fornax. Limits on the number of 7-ray events above the minimal energy 
Ejnin have also been computed (see Table H) and are used in Section |4| for the calculation 
of upper limits on the 7-ray flux. 



4. 7-ray flux upper limits 



Upper limits on the number of observed 7-rays above a minimal energy Emm can 
be translated into an upper limit on the observed 7-ray flux O-y if the energy spectrum 
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^max 


NoN 


Nqff 




Significance 


0.1° 


160 


122 


71 


2.3 


0.5° 


3062 


2971 


243 


1.2 


1.0° 


11677 


11588 


388 


0.6 



Table 3: Numbers of VHE 7-ray events from the direction of the Fornax galaxy cluster cen- 
tre, using three different opening angles for the observation. Column 1 gives the opening 
angle 6max, columns 2 and 3 the numbers of 7-ray candidates in the ON region, Nqn, and 
the normalized number of 7-ray in the OFF region, Nqpf/ respectiv ely. Column 4 gives the 
95%C .L. upper limit on the number of 7-ray events according to iFeldman and Cousins 
(|l998h . The significanc e of the nu r nbers of 7-ray candidates in the ON region is stated in 
column 5 according to lLi and Mai (|l983h . 



dNy/ dEy of the source is assumed to be known, as indicated by equation |9l 



dE 



dNy , 

'^dl;^^^^ 



^obs 



dEj Aeff(E^)-^^(E^^ 



(9) 



Here, Tobs and Agff denote the target observation time and the instrument's effective 
collectio n area, respectively. The in trinsic spectra of standard astrophysical VHE 7-ray 
sources ( Hinton and Hofmann 2009}) typically follow power-law behavior of index T ~ 
2 — 3. Upper limits at 95% C.L. on the integral flux above the minimum energy (cf. Sec- 
tion [S]) are given in Table |4] for different source spectrum indices. 

Dark matter annihilation spectra depends on the assumed annihilation final stat es of 



the D M model. For instance, some supersymmetric extensions of the Standard Model (|Tungman et al. 
19961) predict the neutralino as the lightest stable supersymetric particle, which would be 
a good dark matter candidate. In general, the self-annihilation of neutralinos will give 
rise to a continuous 7-ray spectrum from the decay of neutral pions, which are pro- 
duced in the hadronisation process of final-state quarks and gauge bosons. Universal 
extra-dimensional (UED) extensions of the SM also provide suitable DM candidates. In 
these models, the first Kaluza-Klein (KK) mode of the hypercharge ga uge boson B^^) is 
the l ightest KK particle (LKP) and it can be a DM particle candidate (|Servant and Tait 
20031) . Nevertheless, in the absence of a preferred DM particle model, constraints are 
presented here in a model-independent way, i.e. for given pure pair annihilation final 
state for the DM pair annihilation processes and DM particle mass. The only specific 
DM particle model studied here is the KK B^^^ particle model, where the branching ra- 
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tios of each annihilation channel are known. A wide range of dark matter masses is 
investigated from about 100 GeV up to 100 TeV. A model-independent upper bound 
on the dark matter mass car \ be derived from unitarity for t hermally produced DM as 
done in the semina l pap er of iGriest and Kamionkowskil (|1990( ) and subsequent studies of 
Beacom et al. ( 2007) and Mack al. ( 20081) . Assuming the current DM relic density mea- 
sured by WMAP ([Larson et fl/.ll201ll) , the inferred value is about 100 TeV. Figure |3] shows 
different annihilation spectra for 1 TeV mass dark matter particles. Spectra of DM par - 
ticles annihilating i nto hh, W"'"W~ and r^r ~ pairs are extracted from ICirelli et al\ (|201l[ ), 
and calculated from IServant and TaitI ( 2003 ) for Kaluza-Klein B^'^^ annihilation. Flux up- 
per limits as function of the DM particle mass are presented in Figure |4] assuming DM 
annihilation purely into bb , W~^W~ and t~^t~ and an opening angle of the integration of 
0.1°. Flux upper limits reaches 10 



"12 cm-2 s-i for 1 TeV DM mass. 



^max 


jy95/oC.L.^£^ > Emin) 


r = 1.5 


>Einin)(10-12 cm-2s-l) 

T = 2.5 


0.1° 


41.3 


0.8 


1.0 


0.5° 


135.1 


2.3 


3.3 


1.0° 


403.5 


6.8 


10.0 



Table 4: Upper limits on the VHE 7-ray flux from the direction of Fornax, assuming a 
power-law spectrum with spectral index T between 1.5 and 2.5. Column 1 gives the 
opening angle of the integration region 6max, column 2 the upper limits on the number 
of observed 7-rays above the minimum energy Emin = 260 GeV, calculated at 95% C.L.. 
Columns 3 and 4 list the 95% C.L. integrated flux limits above the minimum energy, for 
two power law indices. 



Recent studies (leltema et flZ.ll2009l:IPinzke et flZ.l2(ffl9l:IPinzke and Pfrommerlboioh have 
computed the cosmic-ray induced 7-ray flux from pion decay s using a cosmological sim- 
ulation of a sample of 14 galaxy clusters (jPfrommer et al. 1120081) . Since the electron induced 
7-ray flux from inverse Co mpton is found to b e systematically subdominant compared to 
the pion dec ay 7-ray flux (|][eltema et fl/.ll2009r) . this contribution is not considered. Using 
the results of Pinzke et al. (I2OO9I) . the 7-ray flux above 260 GeV for Fornax is expected to lie 
between a few 10"^^ cm"^ s"^ and 10"^^ cm"^ s"^ for an opening angle of observation of 
1.0°. The flux is about 2-to-3 orders of magnitude lower than the upper limits presented 
in Table m thus this scenario cannot be constrained. 

Assuming a typical value of the annihilation cross-section for thermally-produced 
DM, {av) = 3 X lO'^*" cm'^s"^, a mass of 1 TeV and the NFW profile of DM density profile 
of Fornax RB02, the predicted DM 7-ray flux is found to be a few 10"^'^ cm"^ s"^. This 
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estimate takes into account the 7-ray enhancement due to dark halo substructure and 
the Sommerfeld enhancement (see section [S]) to the overall DM 7-ray flux. Therefore the 
dominant 7-ray signal is expected to originate from DM annihilations. Constraints on the 
DM-only scenario are derived in the following section. 



5. Exclusion limits on dark matter annihilations 

Upper limits at 95% C.L. on the dark matter velocity-weighted annihilation cross- 
section can be derived from the following formula: 

^ ^ T„bs/(An)An dN^(E^) • ^ 



— r-eiiv-r. ^£ 



7 



The factor / is extracted from Section |2l The exclusion limits as a function of the DM 
particle mass ^dm for different DM halo profile models are depicted in Figures |5] and [6] 
for DM particles annihilating exclusively into bb and B^^^ particles, respectively. Predic- 
tions for (av) as function of th e B*-^-* particle mass are given in Figure [6] within the UED 



framework of IServant and TaitI (|2003h . As an illustration of a p ossible change in this p re 



diction, a range of predicted {(Tv) is extracted from Figure 2 of lArrenberg et al.\ (j2008|) , in 
the case of a mass splitting between the LKP and the next lightest KK particle down to 1%. 
In the TeV range the 95% C.L. upper limit on the annihilation cross-section (av) reaches 
10~^^ cm^s~^. Exclusion limits as a function of the DM particle mass mdm assuming 
DM particle annihilating into bb, t~^t~ and W^W~ are presented in Figure[Zlfor the RB02 
NIW profile. Stronger constraints are obtained for masses below 1 TeV in the t~^t~ where 
the 95% C.L. upper limit on (av) reaches 10~^^ cm'^s"^. The Fermi-LAI e xclusion limit for 



Forna x is added in Figure |5] (pink dashed-line), extending up to 1 TeV (jAckermann et al. 



2010a[ ). It is based on the RB02 NEW profile and a 7-ray spectrum which assumes anni- 
hilation to bb pairs. Below 1 TeV, the Fermi-LAJ results provide stronger limits than the 
H.E.S.S. results. However, the H.E.S.S. limits well complement the DM constraints in the 
TeV range. 

Other DM particle models give rise to modifications of the 7-ray annihilation spec- 
trum which may increase the predicted 7-ray flux. Some of them are considered in the 
following. 
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5.1. Radiative correction: Internal bremsstrahlung 



In the annihilation of dark matter particles to charged final states, internal bre msstrahlung 
processes can contribute sig nificantly to the high-energy end of the 7-ray spectrum (jBergstrom et al. 
2OO5I: iBringmann et fl/.ll2008|) . Adding this effect to the continuous spectrum of secondary 
7-rays from pion decay, the total spectrum is given by 



7^7 



dE 



7 



dE 



(11) 



7 



The r nagnitude of this effect depends on the intrinsic properties of the dark matter par- 
ticle. IBringmann et al\ (l2008h pro vide an approximation that is valid for wino-like neu- 
tralinos ( Moroi and Randall boool) . The annihilation spectrum for a 1 TeV wino is shown 
in Figure lU This parametrization is used in the calculation of the 95% C.L. upper limit 
on the velocity-weighted annihilation cross-section as a function of the DM particle mass, 
presented in Figures [71 and HI The internal bremsstrahlung affects the exclusion limits 
mostly in the low mass DM particle regime, where its contribution to the total number of 
7-rays in the H.E.S.S. acceptance is largest. 



5.2. Leptophilic models 



Recent r neasurements of c osmic ele ctron and positron spect ra by PAMELA (lAdriani et al 



2OO9I) ATIC jchang et altoO^ . H.E.S.S. (lAharonian et al. l2009bl) and Fermi-LAT (lAckermann et al.\ 



2010cl) have been explained in terms of DM annihilation primarily into leptonic final states 



(to avoid an over-pro duction of anti-protons), hereafter referred to as leptophilic models. 
Bergstrom et al. show that the Fermi-LAI electron spectrum and the PAMELA ex- 

cess in positron data can be well explained by annihilation purely into fi'^fi" pairs. In this 
scenario, 7-rays are expected from final state radiation (FS R) of the pair. While 

this final state is rarely found in supersymmetric models ( Tungman et al. 19961) , some 
particl e physics models predict the annihilation to occur p redominantly to lepton final 
states (|Arkani-Hamed et al\\200% iNomura and Thaler! 120091) . The subsequent muon de- 
cay into positrons and electrons may lead to an additional 7-ray emission component by 
Inverse Compton (IC) up-scattering of background photons, such as those of the cosmic 
microwave background (CMB). If the electron /positron energy loss time scale is much 
shorter than the spatial diffusion time scale, the IC contribution to the 7-ray flux may 
be significant. In galaxy clusters, the energy loss term is dominated by the IC compo- 
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nent (|Colafrancesco et alluOOw . The total 7-ray spectrum is then given by 



7^7 



dE 



7 



dE 



(12) 



7 



After extracting the FSR parametrization from lBovyl (|2009l) , the IC c omponent of the anni 
hilation spectrum was calculated following the method described in lProfumo and Teltema 
(|2009|) . The total annihilation spectrum for a 1 TeV dark matter particle annihilating to 
pairs is shown in Figured The energy of the IC emission peak is driven by 
electrons /positrons of energy Ee ^ mDM/2 up-scattering target photon s in a radiation 
field of average energy e = 2.73 K and is given by E}^ ^ e(Ee/me)^ (|Longairl ^992^ . 
Consequently, the enhancement of the 7-ray flux in the H.E.S.S. energy range is found to 
lower the exclusion limits only for very high DM masses, > 10 TeV. The limits are 

enhanced by a factor of ^10. The Fe rmi-LAJ exclusion limit for Fornax is added (gray 
dashed-line), extending up to 10 TeV (jAckermann et flZ.ll2010al ). Due to the IC component, 
below a few tens of TeV the Fermi-LAT results provide stronger limits than the H.E.S.S. re- 
sults. However, since for DM particle masses above 10 TeV the IC emission peak falls out 
of the Fermi-LAT energy acceptance, the IC spectra becomes harder in the same energy 
range. The Fermi-LAT limits for DM particle masses above 10 TeV would tend to raise 
with a stronger slope than the slope in between 1 and 10 TeV. Thus H.E.S.S. limits would 
well-complement the Fermi-LAT constraints in the DM mass range higher than 10 TeV. 

7-rays from IC emission are also expected in the case of DM particles annihilating 
purely into bb. In the H.E.S.S. energy range for high DM masses (> 10 TeV) annihilating 
in the bb channel, the expected num ber of 7-rays inclu ding IC emission is lower than in 
the fi~^pi~ channel (see, for instance, ICirelli et alWlOUi) . This qualitative estimate in the 
Fermi-LAT energy range (80 MeV - 300 GeV) shows that the number of expected 7-rays 
including IC emission for DM particle masses between 1 and 10 TeV is lower in the bb 
than in the channel by at least a factor of 2. Since the (av) exclusion limits are 

roughly scaled by the number of expected 7-rays , a qualitative estimate of the Fermi-LAT 
limits including the IC component in the bb channel should not be better than their limits 
in the channel. 



5.3. Sotnmerfeld enhancement 



The self-annihilation cross-section of dark matter particles can be enhance d with re- 



spect to its value (crv )^ during thermal freeze-out by the Sommerfeld effect (see e.g. lHisano et al. 



2004: |Profumdl2005|) . This is a velocity-dependent quantum mechanical effect: If the rel- 



ative velocity of two annihilating particles is sufficiently low, the effective annihilation 
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cross-section can be boosted by multiple exchange of the force carrier bosons. This can be 
parametrized by a boost factor S, as defined by: 



{(TV 



eff 



S X (crv) 



(13) 



Lattanzi and Silkl (|2009l) consider the case of a Sommerfeld boost due to the weak force 
which can arise if the dark matter particle is a wtno-like neutraltno. As a result of the 
masses and couplings of the weak gauge bosons, the boost is strongest for a DM particle 
mass of about 4.5 TeV, with resonance-like features appearing for higher masses. This 
effect was proposed to account for the PAMELA /ATIC data e xcess, where a boo st of 10^ 
or more is required for neutralinos with masses of 1-10 TeV (jCirelli et alluOO^ . It was 
shown that the b oost would be maximal in the dwarf galaxies and in their substructures 
(jPieri et flZ.ll2009h , due to the low DM particle velocity dispersion in these objects. 



In the Fornax galaxy cluster, the velocity dispersion and hence the mean relative 
velocity of "test masses" such as star s, globular clusters or galaxies is of the order of a 
few 100 km s"^ jSchuberth et fl/.lboioL hence j6 = ^ 10~^. Assuming that the same 

velocity distribution holds true for DM particles, limits on [av) S were derived which 
are shown in Figure [8] for a signal integration radius of 1.0° and the RB02 NFW profile. 
Although the DM velocity dispersion is about one order of magnitude higher than in 
dwarf galaxies, a boost of ^10^ is obtained for DM particle masses around 4.5 TeV. The 
resonance-like feature is clearly visible for masses above 4.5 TeV. Outside the resonances, 
the limits on (av)^^^/ S are tightened by more than one order of magnitude for dark matter 
particles heavier than about 3 TeV. 



5.4, Enhancement from dark matter substructures 

The effect of DM substructures inside the opening angle of 0.1° and 1.0° are presented 
in Figure HI using the enhancement values calculated in Section 12.21 The enhancements 
to the 95% C.L. upper limits on (av) are estimated using the two limiting masses of sub- 
structures Miijn. In the TeV range, the upper limit on (av) is at the 10~^^ cm-^s"^ level. 
The joint enhancement due to the Sommerfeld effect added to the IB and the substruc- 
tures contribution is plotted in Figure HI In the most optimistic model, with the largest 
enhancement by substructures and the Sommerfeld effect, the 95% C.L. upper limit on 
((rz;)gff reaches 10~^^ cm^s~^, thus probing natural values for thermally-produced DM. 
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6. Summary 

The Fornax galaxy cluster was observed with the H.E.S.S. telescope array to search 
for VHE 7-rays from dark matter self-annihilation. No significant 7-ray signal was found 
and upper limits on the 7-ray flux were derived for power-law and DM spectra, at the 
level of 10"12 cm-^s-i above 260 GeV . 

Assuming several different models of particle dark matter and using published mod- 
els of the dark matter density distribution in the halo, exclusion limits on the DM self- 
annihilation cross-section as a function of the DM particle mass were derived. Particular 
consideration was given to possible enhancements of the expected 7-ray flux which could 
be caused by DM halo substructure or the Sommerfeld effect. For a DM mass of 1 TeV, 
the exclusion limits reach values of {av) ~ 10~^^ — 10~^^ cm'^s"^, depending on DM 
model and halo properties, without the substructures contribution, and {av) ^ 10"^"^ — 
10 cm^s~^ when considering the substructures signal enhancement. At Mdm ~ 4.5 
TeV, a possible Sommerfeld resonance could lower the upper limit to 10~^^ cm'^s"^. 



Compared to observati ons of dwarf spheroidal g alaxies (see for instance I Abramowski et at 



2011al) or globular clusters (jAbramowski et alxdQllm , these limits reach roughly the same 



order of magnitude. The choice of different tracers to derive the DM halo profile in Fornax 
galaxy cluster allows to well constraint the uncertainty in the expected signal. The poorly 
constrained, but plausibly stronger subhalo enhancement in the Fornax cluster induces a 
uncertainty in the expected signal of about two orders of magnitude. 

With an optimistic joint 7-ray signal enhancement by halo substructures and the 
Sommerfeld effect, the limits on (crv) reach the values predicted for thermal relic dark 
matter. Additionally, they extend the exclusions calculated from Fermz-LAT observations 
of galaxy clusters to higher DM particle masses. 
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Fig. 1. — Substructure 7-ray flux enhancement as function of the opening angle of inte- 
gration. Two values of the limiting mass of substructures are used: Mii^ = 10~^Mq, for 
the high (HIGH) boost (solid line), and = 5 x W~^Mq, for the medium (MED) boost 
(dashed line). The RB02 profile is chosen as the smooth host DM halo. 
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Fig. 2. — Left: Significance map in equatorial coordinates, calculated according to the Li 
& Ma method (|Li and Mai 1 19831) , with an oversampling radius of 0.1° . The white circle 
denotes the 0.1° integration region. No significant excess is seen at the target position. 
Right: Distribution of the significance. The solid line is a Gaussian fitted to the data. The 
significance distribution is well described by a normal distribution. 
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Fig. 3. — Photon spectra for 1 TeV dark matter particles self-annihilating in different chan- 
nels. Spectra from DM annihilating purely into bb (dot-dashed line), t+t" (black solid 
line) and W^W~ (long-dashed dotted line) are shown. The latter shows the effect of Inter- 
nal Bremsstrahlung (IB) occuring for the W^W~ channel. The 7-ray spectrum from the 
annihilation of B*^^^ hypergauge boson pairs arising in Kaluza-Klein (KK) models with 
UED is also plotted (dotted line). The long dashed line show the photon spectra from 
final-state radiation (FSR) and the inverse Compton (IC) scattering contribution in the 
case of DM particles annihilating into muon pairs. 
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Fig. 4. — Upper limits 95% C.L. on the 7-ray flux as a function of the DM particle mass 
for Ejnin= 260 GeV from the direction of Fornax. DM particles annihilating into bb (solid 
line) , PV+PV" (dotted line) and t~^t~ (dashed line) pairs are considered. 
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Fig. 5. — Upper limit at 95% C.L. on the velocity-weighted annihilation cross-section (av) 
as a function of the DM particle mass, considering DM particles annihilating purely into 
bb pairs. The limits are given for an integration angle 9max = 0.1°. Various DM halo pro- 
files are considered: NFW profiles, SRIO aio (blue solid line), DWOl (black solid line), 
RB02 (pink solid line) and RS08 (green solid line), and Burkert profiles, SRIO a(, (red dot- 
ted lin e) and am (blue solid l ine). See Table [1] for more details. The Fermi-LAH upper 
limits (jAckermann et alWlOlOm for the NFW profile RB02 are also plotted. 
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Fig. 6. — Kaluza-Klein hypergauge boson B^^) dark matter: Upper limit at 95% C.L. on 
{(Tv) as function of the B'^^^ mass towards Fornax. The limits are given for an integration 
angle 0max = 0.1°. The NFW profiles, SRIO aio (blue solid line), DWOl (black solid line), 
RB02 (pink solid line) and RS08 (green solid line), and Burkert profiles, SRIO a6 (red dot- 
ted line) and aio (blue solid line). See Tabled] for more details. The prediction of {av) as 
function of the B^^^ mass is given (dotted-line). A range for this predictions is given in 
case of a mass splitting between the LKP and the next LKP down to 1% (dashed area). 
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Fig. 7. — The effect of different DM particle models: Upper limit at 95% C.L. on (crv) as 
function of the DM particle mass. The limits are given for ^max = 0.1°and the NFW profile 
RB02. The limits are shown for DM particles annihilating into hi (gray solid line) , W"*" W~ 
(gray dash-dotted line), t'*"t~ (gray long-dash-dotted line) pairs. The effect of Internal 
Bremsstrahlung (IB) occuring for the W+W~ channel is plotted in gray long-dashed line. 
The black solid line shows the limits for DM annihilating into p airs including the 



effect of inverse Compton (IC) scattering. The Fermi-LAT upper limits (|Ackermann et al 
2010al) for the NFW profile RB02 and for an DM annihilating into J4~^ji~ pairs including 



the effect of IC scattering are also plotted (black dotted line). See section 12.21 for more 
details. 
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mDM(TeV) 

Fig. 8. — The Sommerfeld effect: Upper limits at 95% C.L. on the effective annihilation 
cross-section {o'v)^^^ = {crv)Q/S as a function of the DM particle mass annihilating into 
W pairs. The black line denotes the cross-section limit for ^max = 1-0° without 7-ray flux 
enhancement, the dashed blue line shows the effect of halo substructure (using the "high 
boost", cf. Fig. H). The solid green and blue lines show the limit for the case of Wino 
dark matter annihilation enhanced by the Sommerfeld effect, with and without including 
Internal Bremsstrahlung, respectively. The DM halo model RB02 is used (see Tableland 
main text for more details). A typical value of the annihilation cross-section for thermally- 
produced DM is also plotted. 
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10' 1 10 10^ 

mDM(TeV) 

Fig. 9. — The effect of DM halo substructures: Upper limit at 95% C.L. on {av) as function 
of the DM particle mass annihilating purely into hh pairs. The limits are given for 0max = 
0.1° (dashed lines) and 0max = 1-0° (solid lines). The DM halo model RB02 is used (see 
Table [Hand main text for more details). In addition, the effect of halo substructures on the 
{(Tv) limits is plotted. The "medium boost" (MED) with = 5 x 10"^ M,, (blue lines) 
and the "high boost" (HIGH) with Mum = 10~^ Mq (red lines) are considered. 
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